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The non-equilibrium quantum field dynamics is usually described in the closed-time-path for- 
malism. The initial state correlations are introduced into the generating functional by non-local 
source terms. We propose a functional approach to the Dyson-Schwinger equation, which treats 
the non-local and local source terms in the same way. In this approach, the generating functional 
fNj , is formulated for the connected Green functions and one-particle-irreducible vertices. The great 

advantages of our approach over the widely used two-particle- irreducible method are that it is much 
simpler and that it is easy to implement the procedure in a computer program to automatically gen- 
s ! | erate the Feynman diagrams for a given process. The method is then applied to a pure gluon plasma 

to derive the gauge- covariant transport equation from the Dyson-Schwinger equation in the back- 
ground covariant gauge. We discuss the structure of the kinetic equation and show its relationship 
with the classical one. We derive the gauge-covariant collision part and present an approximation 
in the vicinity of equilibrium. The role of the non-local source kernel in the non-equilibrium system 
is discussed in the context of a free scalar field. 



PACS numbers: 12.38.Mh, 25.75.-q, 24.85.+p, 11.15.Kc 
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I. INTRODUCTION 

o 

The ultimate goal of ultra-relativistic nucleus-nucleus collisions is to study the properties of strongly interacting 
matter under extreme conditions of high energy density Qj. Quantum Chromodynamics (QCD) predicts that strongly 
interacting matter undergoes a phase transition from a state of hadronic constituents to a plasma of unbounded quarks 
O • and gluons (QGP) [pf. The QGP is considered as a partonic system being at (or close to) local thermal equilibrium. 
Thus, to study the conditions for the possible formation of a QGP in heavy-ion collisions one needs to address the 
question of thermalization of the initially produced partonic medium || ||. There are various theoretical predictions 
suggesting that the appearance of the QGP could modify the properties of physical quantities that are measured in 
heavy-ion collisions. The suppression of charmonium production pi was argued to be a possible consequence of the 
collective effects in the thermalized and deconfined medium. The jet quenching [7lwas predicted to be due to 
the radiative energy loss of partons penetrating the QGP. Strangeness [p and dilcpton |p production yields could be 
$_i modified in the thermal QCD medium. A variety of models for the initial conditions in ultra-relativistic heavy-ion 
collisions suggest that at the early stage the medium is dominated by the gluon degree of freedom 10 . The transport 
equation for a pure gluon plasma is thus of special interest. 

The usual treatment of the gluon transport equation is based on the decomposition of the gluon field into a mean 
field and a quantum fluctuation. Under this approximation the gluon transport equation then describes the kinetics 
of the quanta in the classical mean field jn], [l2[ This picture is somewhat similar to what was used while 

studying the energy loss of the fast parton moving in the soft mean field || f?J. To include the classical chromofield 
into QCD in a proper way, one uses the background field method of QCD (BG-QCD) introduced by DeWitt and 't 
Hooft jli], [ll| |1| [It], [IS). The advantage of BG-QCD is that it is formulated in an explicit gauge-invariant manner. 
The BG-QCD is a very suitable method to describe the properties of a QCD medium created in the initial state of 
heavy-ion collisions. The time evolution of a quantum system being off equilibrium can be, in principle, obtained by 
solving the Dyson-Schwinger equation (DSE) denned on a closed-time-path (or the Kadanoff-Baym equations). If the 
kinetic scale, describing long-range correlations, is much larger than the scale of quantum fluctuations, the DSE may 
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be reduced into a much simpler form of the transport equation by a gradient expansion [[L9| |2(], |2l|, [22j . To derive the 
transport equation in the presence of a classical chromofield in an explicit gauge-invariant or covariant way, one thus 
combines the BG-QCD method and the closed-time-path (CTP) formalism. 

To our knowledge, the first study of the kinetics of a classical particle with non-Abelian gauge degrees of freedom 
propagating in a non-Abelian classical gauge field is due to Wong p3[ . There were many efforts in the literature to 
derive the transport equation for the QCD medium |2f| ^6j, particularly by using the BG-QCD method [^7], |2^, g9[ . 
To our knowledge, the first application of this method was done by Elze p7[ | to derive the transport equation for gluons. 
His approach was based on the Yang-Mills field equation and the second quantization in the operator representation 
of the quantum field theory. However the equation obtained is in a complicated form. In Q , the transport equation 
for the gluon was derived in the conventional QCD by using the light-cone gauge and CTP formalism. The gluon field 
was decomposed into a hard and a soft part treated as the classical field. This decomposition, however, was not done 
in a gauge-covariant way. Most recently Blaizot and Iancu have derived the Boltzmann equation for the QCD plasma 
p8[ , applying the CTP formalism and the BG-QCD method which guarantees the gauge-covariant decomposition of 
the gauge field. 

To derive the Boltzmann equation one needs, in general, to make a set of approximations. This usually involves 
a gradient expansion of the DSE and the perturbative derivation of the collision term. In Ref. p8| an additional 
assumption has been made that the system is close to equilibrium. Consequently, the transport equation for the QCD 
plasma obtained in Ref. |28| was linearized with respect to the off-equilibrium fluctuations. 

In this paper we propose a different approach to derive the Boltzmann equation for the gluon plasma in the BG-QCD. 
First we apply the functional approach to the DSE that treats the non-local and local source terms in the same way. 
We strictly stick to the functional definition of the one-particle-irreducible (1-PI) vertices and the connected Green 
functions (CGF). Furthermore, we use the DeWitt notation, which, in our opinion, results in a simple structure for 
the generating functionals. The current approach has a great advantage over the widely used two-particle-irreducible 
method that it is much simpler and can automatically generate all necessary Feynman diagrams. 

In a heuristic discussion on the role of the non-local kernel for a free scalar field, we show that if the initial time 
is in the remote past the kernel provides only a correction to the homogeneous solution of the DSE and preserves 
its structure. We also see in this simple model that if the initial time is not in the remote past, the non-local kernel 
brings the time dependence to the DSE and breaks the assumed structure of the homogeneous solution. 

From the DSE, we derive the transport equation in a gauge-covariant way. Our derivation is quite general as it 
does not require any additional assumptions, such as a special form for the gauge-covariant Green function (GF) or 
that the system is near equilibrium. Consequently our equation is not linearized with respect to the off-equilibrium 
fluctuations. We use the background covariant instead of the background Coulomb gauge as was used in Ref. [ p8| . 
Therefore our results preserve an explicit Lorentz covariance and have a compact structure. However, we have to 
include the ghost fields to cancel the non-physical degrees of freedom of the gauge field. We note that the resulting 
kinetic equation has a structure similar to the one previously derived in [ f36| , based on p"2| ] in QCD and assuming the 
2-point gauge-covariant GF (or the Wigner function) to be proportional to the quadratic product of the generators 
for the fundamental color representation. 

We discuss the structure of the kinetic part of the transport equation derived here and compare it with the classical 
kinetic equation. In the quantum case the kinetic equation describes the time evolution of the gauge-covariant Wigner 
function, which is a matrix in the adjoint color space. Therefore, it contains many non-Abelian features, which are 
absent from the well known classical equation. However, a notable result is that, as in the classical case, it contains 
a term that corresponds to the color precession |30|. This is the non-Abelian analogue to the Larmor precession for 
particles with magnetic moments in an external magnetic field. We argue that this term is necessary to preserve the 
gauge covariance of the resulting transport equation. We also discuss the structure of the kinetic part if the system 
exhibits only a small deviation from equilibrium. 

We derive the gauge-covariant collision part of the Boltzmann equation and present its linearized form with respect 
to the off-equilibrium fluctuations. Finally applying the transversality condition for the gauge-covariant GF, together 
with some other approximations, the collision part is further simplified and shows the explicit collision and damping 
terms. 

The paper is organized as follows. In Section II we introduce the basic concept of BG-QCD. In Section III we 
present two equivalent, in a path-integral sense, methods to derive the classical equation of motion for the gluon. In 
Section IV, we describe the functional approach to the DSE in the vacuum. We derive the DSE for the 2- and 3-point 
GF in the background covariant gauge. In Section V we present a derivation of the DSE for the non-equilibrium 
gluon plasma within the CTP formalism. Section VI contains a heuristic discussion on the role of the non-local source 
kernel. In Section VII we derive the transport equation from the DSE, applying the gradient expansion. In Section 
VIII the gauge-covariant collision part is derived and its structure under some approximations is discussed. 

In the paper we use — diag(l, —1, —1, —1) as the metric tensor. The Lorentz indices are written as subscripts 
and color indices as superscripts to the relevant quantities. For a pure gluon plasma the color field transforms only 
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under the SU (3) c adjoint representation, thus all color indices are adjoint ones. We use the following notation for the 
gauge field and its field strength tensor: = A^T a and F^[A] = F^[A]T a , where (T a ) i: > = if a: > are the generators 

of the SU(3) C adjoint representation and f abc are the 5t/(3) c structure constants. The two-point GF or self-energy 
(SE) are treated as matrices, thus the color and/or Lorentz indices are sometimes omitted. 

For convenience, we list all abbreviations that we use in this paper: background gauge QCD (BG-QCD), Dyson- 
Schwinger equation (DSE), closed-time-path (CTP), Green function (GF), connected Green function (CGF), one- 
particle- irreducible (1-PI) and self-energy (SE). 



II. BACKGROUND FIELD METHOD 



Any physical quantity calculated in QCD is gauge invariant and independent of the particular gauge chosen. The 
classical Lagrangian for the gauge field exhibits an explicit gauge invariance. However, when quantizing the theory 
in a particular gauge, one introduces the gauge fixing and ghost terms that break the explicit gauge invariance of the 
Lagrangian. The background field method is such a technique that allows a gauge (background gauge) to be fixed 
without losing the classical gauge invariance. 

In BG-QCD the conventional gluon field is expressed as a sum of a classical background field A and a quantum 



fluctuation Q. The action is given by |16, 17, [L8| 

S = So + Sfi x + Sghost + S 

src-i 

S = -~ J d*x(Fi v [A] + ir; AQi - D^AjQf, + gf^Q^Qt) 2 , 
S gh ost = J d A xC l D^[A]D^ k [A + Q]C\ 

Ssrc = J d 4 x(JlQ^ + fc l + CV) . (1) 
The generating functional for the GF reads 

Z[A, J,£,£] = J [dQ][dC][dC] exp(iS), (2) 

where the quantum fluctuations of the gluon field are the integration variables in the functional integral; C l ', C % 

are the ghost and antighost fields; £ , £ l are the external sources coupling to the ghost and antighost fields respectively. 
The covariant derivative D*?'[A] is defined by D^[A] = d^8 lj - ig(T a ) ij 'A*. The field strength tensor F'^[A] for the 
background field is given by F^[A} = d^A\, - d v A\ + gfi k A^A k v . 

There are two types of gauge transformations that leave Z invariant 18 . The type I transformations are given by 



A' = UA^U- 1 + -Ud^U-\ 
9 

Q'^ = uq^u-\ j; - uj^u-\ 
= ucu-\ t = uzu-\ 

C' = UClf- 1 , f = U^U- 1 , (3) 
where U(x) = cxp(iguj a (x)T a ) . For an infinitesimal gauge transformation the fields transform as 





= D%[A]^, 








= u.r jh Qi~ L - 




= gf Jk J^ 


8C l 


= gf ijk C j u k , 


sf = 




5C l 


= gf ijk C j ui k , 


sc = 


gf ijk Z j cj k . 



(4) 

In the type II transformations, the background field does not change, but the gluon field transforms as Q'^ = U(A^ + 
Qf^U^ 1 + ig~ 1 Ud f2 U~ 1 . In this paper, as in Ref. [^7), type I transformations are relevant where the background field 
transforms in a conventional way and the gluon transforms like a matter field. 
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In BG-QCD, we define the generating functional W[A,J,f,f] for the CGF and T[A, (Q), (C), (C)} for the 1-PI 
vertex. Both W and T are gauge-invariant functionals, thus should be invariant under the gauge transformations (|^) 
and (0). Consequently, one can formulate for T/W a Ward identity that corresponds to this invariance Q. 



III. EQUATION OF MOTION FOR GLUON AND CLASSICAL FIELD 



The equation of motion for the gluon field Q is obtained by using the conventional variation approach. Shifting the 
gluon field Q^(x) by a small amount Q l fl (x) = Q^(x) + v^x), the variation of the functional W is 

5W = -idlnZ = Z- 1 J [dQ][dC][dC] J d 4 x-^-v^(x) exp{iS) . (5) 

Because we only shift a variable of integration in the functional integral, there should be no change for W. This 
requirement leads to the equation of motion for Q % - 



SS 



SQl(x) 



Z- 1 J [dQ}[dC}[dC}-^-eM^S) = , (6) 



where the expectation value (• • •) is taken in the path integral sense. In the absence of the external sources J, £ and 
£ the equation of motion becomes 

D i J[A]F^[A]=(j^ i ) , (7) 

with the induced current j is given by: 

j = JO + jgf , 

ft* = -gf abc Dl a [A](Q^ b Q^ c ) - gf lda Qt{D^ ab [A}Q^ b - D^ ab [A]Q^ b ) 
_ g 2pdajabcQdQv,bQn,c _ [A](D v ' jk [A}Q^' k ~ D^ k [A]Q v - k ) 
-gf ijk QlF^ k [A] , 

fj = g f ka C 3 D^ kj \A]C l ~ -D^\A}D^ k [A}Q k (8) 
gj a 

where jo comes from So and j g f from the ghost and gauge-fixing terms. One can verify that j transforms like a matter 
field. 

In fact, we have (5S/SQ) — (5S/SA) = 0. The proof can be found in Refs. 18 and From the condition that 
(5S/5A) = 0, and after setting all external sources to zero, we obtain 

D iJ[A]F^[A} = {j'^) , (9) 

where j' is defined by 



j = JO + jgf , 

J gf 

+-gf l3k Q^ k Dl l [A]Q vl . (10) 



,-'/<■< _ - g f kil C j D "' ki [A]C l - gf ljk & 'D"' kl [A + Q]C l 



Note that j' is different from j because j g f ^ j' However, the expectation values of these two currents are equal. 
Following from the identity [A]Dl k [A]F v ^' k [A] = 0, we also have D^[A](j^) = D l J[A](f^) = 0. 



IV. A FUNCTIONAL APPROACH TO THE DYSON-SCHWINGER EQUATION IN BG-QCD 

In the previous section, we have discussed the classical equation of motion for the background and quantum fields. 
In the following we will show how to obtain the DSE. We will use a convenient functional approach, which has the 
advantage that the non-local and the local terms are treated in the same way. In our derivation we always use the 
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notation introduced by DeWitt |32|, which gives each formula a simple face and a clear structure. In DeWitt notation 
the classical action (|l|) can be written in the following form: 

S = So 4- S src , 

So = -T^ n (A 2 )A m A n + —T^ np (A 3 )A m A n A p + — T < ^ l l Lpq (A 4 )A rn A n A p A q 

^"2^mn(Q )QmQn + Q^iJnpiQ )QmQnQp + ^^inlipqiQ )QmQnQpQq 

+T$ n (AQ)A m Q n + \v^ np {AQ 2 )A m QnQp + \v^ np {A 2 Q)A m A n Q p 



Q^rnnpq(A Q)A m A n A p Q q + Y^ npq (AQ )A m Q n Q p Q q 

-\r^ p q(A 2 Q 2 )A m A n QpQ q 



^^rnniC^CmCn + Y^ np (CC A)C m C n A p + T^! np (CCQ)C m C n Qp 

+F(®l ipq (CCAQ)C m C n ApQ q + ^^ n p q {CCQ 2 )C m CnQpQq , 

Ssrc — JmQm ~t~ ^m^m ^m^m 3 (H) 

where subscripts m,n,p,q represent all necessary indices, including colour, Lorentz and space-time coordinates. The 
fractions in front of some terms are symmetry factors; Tmn P (AQ 2 ), for example, is the bare vertex attaching one A 
and two Q' s - The explicit forms of these bare vertices are given in Appendix A. We use the convention that the 
repetition of two indices stands for a sum or an integral. Note that the definition of Sq has changed from Eq. ([!]). 

From Section III, one knows that the expectation value of the first derivative of S with respect to Q leads to the 
classical equation of motion. In DeWitt notation this derivative is written as: 

' )S = r^,(Q 2 )Q„, + T^ m (AQ)A n , + T^ n , m (AQ 2 )A m ,Q n , 



SQr 

+ y^m>n>m(A 2 Q)A m >A n i + ^r|J n , p , m (A 3 Q) A m > A n 'Ap> + 7j^ m ) n , plm (A 2 Q 2 )A m >A n >Qpi 



2 m'n'm\ ^ / " fc 6 mn P mv ' 2 

'T^inn'p'iQ )Qn'Qp' + Q^Un'p'q'iQ )Qn'Qp'Qq' + 7^m'n'p'm(AQ )A m 'Qn'Qp' 

,'p 



+T ( % n , m {CCQ)C m ,C n , + T { % n , lm {CCAQ)C m ,C n ,Ap, 



+T'Un'p'm( ( ~'CQ 2 )C m >C n rQp> + J m . (12) 

The classical equation of motion then reads 

SS \ n / SS \ . M „. 

or ( — — ) = - J m . (13) 



SQm / \SQ 

The explicit form of the equation of motion in the absence of J m is given by Eq. (|t]). From the 1-PI generating 
functional, we also have ST/S(Q m ) — —Jmi consequently we have 

SS \ ST 



SQ m I S(Q 

From now on we will always refer to Sq instead of S 1 , we therefore drop the subscript and simply denote it as S. 
Taking the derivative of Eq. (|lj) with respect to (Q p ), we have 

<5 2 r sj n s I ss \ s 2 r s / ss 



S(Q P )5(Q m ) S(Q P ) SJ n \SQ m / S(Q P )S(Q„) SJ n \5Q T 

= -zr pn (g 2 )[(^Q„)-(Q„)a»>] , (15) 

where Z^, is given by Eq. ([l2]), except for the constant terms A, A A, AAA and J, which do not contribute to the 
DSE. The explicit expression for {l} n Q n ) reads: 

(I'M = r^(Q 2 )(Qn>Qn)+r ( °] n , m (AQ 2 )A m ,(Q n ,Q n ) 
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+ 2^m'n'p'm(^ 2( 3 2 )^ m '^ n ' (Qp'Qn) 



2 

+ 2^rn'n'p'm(AQ 3 )A m i (Qn'Qp'Qn) 

+r i ° ) , n , m (CCQ)(C m ,C n ,Q n ) + Y^, n , plm {CCAQ)A pl {C m ,C n ,Q n ) 
+^L P - m (CCQ 2 ){C m ,C n/ Q p ,Q n ) . (16) 

We see that all terms in Eq. (fL5|) are n-point correlation functions, or GFs, and can be expressed in terms of CGFs. 



The relations between the GF and the CGF are given by Eqs. (Bl) and (B2) in Appendix B. A higher-rank CGF can 



be related to lower-rank CGFs and 1-PI vertices by various identities. These identities are summarized in Eqs. (B3) 



and (B4) in Appendix B. For our further discussion we also introduce the following notions and conventions for the 



GF and CGF. We denote the 2-, 3- and 4-point CGF for Q G mn (Q 2 ), G mnp (Q 3 ) and G mnpq (Q 4 ), respectively. For 
convenience, sometimes, we also use a short-hand notation (inn) = G mn (Q 2 ) 1 (mnp) = G mjip (Q 3 ), etc. The 2-, 3- 
and 4-point CGFs for Q and C/C are denoted by G mn (CC), G mnp (CCQ) and G mnpq (CCQ 2 ), respectively. We also 
use the following simplified notations: ([mn]) = G mn {CC), ([mn]p) = G mnp (CCQ), ([mn]pq) = G mnpq (CCQ 2 ) etc. 

We treat A and Q as two independent variables. In the last step of the calculations, we always set (Q) and all 
sources to zero. As it is only in the last step that (Q) is set to zero, the external sources and background field A are 
independent of each other in the intermediate steps. This is the difference between our approach and that used in 
Ref. p7| . Thus, in deriving the DSE for the 2-point GF, we can drop all terms in Eq. ( fl5| ) that are proportional to 
(Q) as we need not take further derivative with respect to (Q). Note that the correlation function (Q m Q n ) becomes 
simply G mn (Q 2 ) after dropping the (Q) term. 
Following the identities 

S 2 W S 2 T 
SJ P SJ q ' 6(Q p )S(Q m ) ~ ~ qm ' 
S 2 W 

JsjpJT, = Gpq ' (17) 



we obtain the DSE for 2-point GF from Eq. (|15|): 

*G4(Q 2 ) = iG( ] mp [Q 2 } + n mp (Q 2 ) , (18) 



where 



5 2 T 

lG "^ (Q2) = 8(Q p )6(Q m ) ' 

n (o 2 ) - n (1) +n (2) 
i^\{Q 2 )G pq = ir^, p ,(g 3 )(Q„,Q p ,g g ) + ir^ Vg ,(Q 4 )(Q I1 ,g p ,Q 9 ,g 9 ) 

+ ^m > 'n'p'm{AQ' i )A m i (Qn'Qp'Qq) j 

n£i(Q 2 )G pq = T% n , m (cCQ)(C m .C nl Q q ) + T ( °l >p> jCCAQ)A p ,(C m ,C n ,Q q ) 

+r i °} n , plm (CCQ 2 )(C ml C n ,Q p/ Q q ) . (19) 

Note that iH pm (Q 2 ) is the normal SE for Q; iH^ is the SE from the gluon loop and ill^ that from the gluon and 
ghost loops. Various bare vertices r(°) can be derived from the classical action ([!]) by taking the functional derivatives 
with respect to the field expectation value. The results are summarized in Appendix A. 



The DSE (18) can be written in an alternative form: 

*^(0)mp^'p9 UmpCpq = lS m q , (20) 
^nm^fQJ m p ~t~ Gjimllmp i^rvp • (21) 
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where we omit the arguments Q 2 in G mp {Q 2 ), G mp (Q 2 ) and iTl pm (Q 2 )- We take Eq. ( |2C| ) as an example to expand 
and analyse IlmpGpq. We know that H m pG pq = TlmpG pq + TlmpG pq where Hmp^G pq are given by Eq. (|l9|). Using 
Eqs. (Bl) and (B3) from Appendix B, we can expand HmpG pq as follows: 

n&Gp, = irW v (Q8 )(n y 3)+ l r W n#j) , m( ^ga )i4ml(n y 9) 

+^ril p , 9 ,(g 4 )[zr Ws , t , u ,(Q 4 )(rV)( s V)(i'g')( w 'g) 

+3zr r ^^,(Q 3 )(r'n' g )( S V)(i' (? ') + 3(nV)(g'g)] . (22) 

Note that the fraction in front of each term is a symmetry factor, which is automatically reproduced. In Fig. 1 we 
show the Feynman diagrams corresponding to the above equation. The solid line there, with two dots at its ends, 
represents a 2-point full GF. The 3-point full GF is drawn as a shaded circle with three solid-line legs, which have 
dots at their outer ends. The 1-PI vertex is a shaded circle with all legs amputated; it has dots on the circle that 
stands for points to which amputated legs are attached. In the same way, we expand Tl mp G pq , which is associated 
with the ghost loop as 



n^G OT = rW n , m (CCQ)([m'n']g) + r^^ p , m (CC^Q)^([m'n']g) 

+^l llp , m {CCQ 2 ){{[ m 'n']p' q ) + ([m'n'Wq)} 

= rW n , ro (^CQ)([mV] ? )+rW nVro (^CAQ)^([m'n'] ? ) 

+r^ n> - m (CCQ 2 ){zr r , s ^(CCQ 2 )([m' S '] g )([r'n'])(tV) 

+iT r , slt ,(CCQ 2 ){m l s')([r'n'] q ){t'p') 
■filWf (CCQ 2 ) ( [m V] ) ( [r V] ) (t'p'q) 
+iT r , s , Vu , (CCQ 2 )( [mis'] ) ( [r'n'} ) (t'p')(u'q) 

+([m'n'])(p'q)} . 



(23) 



In the expansion of Eqs. (|22h and (|23|), we have set (Q), (C) and (C) to zero. Feynman diagrams corresponding to 
Eq. ( p3| ) are shown in Fig. 2. 

Until now we have derived the DSE for the 2-point GF. In the following we extend our analysis to the DSE for the 
3-point GF. Taking the derivative with respect to (Qi) in Eq. (|lH), we obtain 

S 2 T 



-r pn (Q 2 )T u >(Q 2 )((i m QnQi>) - (Ql>)(I m Q: 
+iT lpn (Q 3 )(Q n )(I m )+T pn (Q 2 )T lu (Q 2 ) 

■ \l m )((QnQl') ~ (Qn)(Ql')) + (Qn)((I m Ql>) ~ (OiQl'] 

After setting (Q) = 0, we have 

IWQ 3 ) = -iT lpn (Q 3 )(I m Q n ) -r pn (Q 2 )T w (Q 2 )(I m Q n Qi>) 

+Tp n (Q 2 )T W (Q 2 )(I 1 n [)(QnQl>) , 



(24) 



(25) 



where I m is as I m without the (Q) terms. From Eq. ([l5|), by dropping the term (<5„)(/^ l ), the DSE becomes 
(I m Qn) = ifimn- We apply it to the above equation and obtain 



(I m QnQp) = (QnQp)(I 1 n ' l ) ■ 

Using Eq. (p^|), we explicitly write (I m Q n Q P ) as follows: 

{ilQuQp) = r^l(Q 2 )(Q n ,Q n Q p )+T^ n ,JAQ 2 )A m ,(Q n ,Q n Q p ) 



(26) 
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+ IjF^m' ri' rn(^ Q 2 ) Am' A n ' {Q p i Q n Qp) 

+l^l' P '(Q 3 )(Qn'Q P 'Q n Q P ) + ir^, p , g ,(Q 4 )(Q„,Q p ,Q 9 ,Q n Q p ) 

+ ^ly m W 3 )An' (Qn'Qp'QnQp) 

+rS nVni (^CAQ)^^C n ,Q n Q p ) 

+r£l Vm (CCQ 2 )<C m ,C n ,Q p ,Q„Q p ) . (27) 
We can prove that the disconnected part of (I^QnQp) cancels (Q n Q P ) (!„[), hence we write Eq. ( f26| ) as 

*^(0)mn' [Q 2 ]Cn'np( ( 3 3 ) = — Primp , (28) 

where the 1-PI vertex r mnp is defined by 

Tmnp — 2^ mn 'p'(^ ){Qn'Qp'QnQp) c ~ m'n'p'rn^^^ )A m ' (Qn'Qp'QnQp) c 
"F "O^mn'p'q' ){Qn'Qp'Qq'QnQp) c 



•pi c 



+r!$ n , m {CCQ)(C m ,C n ,Q n l 

~^~^rn'n'p'm (C'CAQ)Ap' (C m 'C n 'Q n Qp) c 

+r n 3 X l ,p, m (CCQ 2 )(C m ,C n/ Qp,Q n Qp) c , (29) 

and the subscript c stands for the connected part. 

We can expand T mnp term by term. The first term |r^, , {Q 3 )(Q n 'Q P 'QnQp)c in Eq. (|||) can be obtained by 



expressing (Qn'Qp'QnQp) by the CGFs of the lower rank. This relation is given by Eq. flB5|) in Appendix B. We can 
identify in this relation disconnected GFs (n'p')(np). After dropping it, we get the connected part (Qn'Qp'QnQp) c- 
The corresponding Feynman diagrams are shown in Fig. 3. The second term in Eq. (|2^) is the same as the first 
one, except that there is an additional A field attached to the bare vertex. The third, fourth and fifth terms can 



be expressed in terms of the lower rank CGFs. The resulting relations are given respectively by Eqs. (B6-B1C) in 
Appendix B. Their corresponding Feynman diagrams are shown in Figs. 4-6. 

In this section, we have derived the DSE for the 2- and 3-point GFs. We strictly stick to the functional definition 
of the 1-PI vertex and the CGF and use the DeWitt notation. Finally, the relations between CGFs and 1-PI vertices 
were recursively applied to express a higher rank CGF in terms of the lower rank ones and 1-PI vertices. The current 
approach has the advantage that the non-local terms can be treated in the same way as the local ones. The difference 
between the local vertex and the non-local one is that the former has a sufficient number of (5-functions to ensure that 
vertex is at the same space-time point. Another advantage of the current approach is that it can produce all needed 
Feynman diagrams automatically. Hence it is easy to implement our approach in a computer algorithm, which can 
automatically generate Feynman diagrams for a given process. 



V. DYSON-SCHWINGER EQUATION IN CLOSED-TIME-PATH FORMALISM 

The non-equilibrium dynamics is usually described in the CTP formalism |3^, |34j . In this section, we will formulate 
the DSE in this formalism. The generating functional Z[J,£,£\ in the CTP formalism reads: 

Z{A ± ,J±,U,I±} = J{dQ + }{dQ_}{dC + }{dC-}{dC + }{dC-} 

■ exp{iS+ - iS- + iK(A± 7 Q±)} 
S+ = S(A + ,Q + ,C + ,C + ,J + ,Z+,Z + ) 

S- = S(A-,Q-,C-,C-,J-,£-,£_) , (30) 

where the classical action S of BG-QCD is given by Eq. (|l|). We denote the total action as Sctp = S*+ — 5_. We 
have omitted the kernel K(A±, Q±). This is because, as we will see in the next section, the kernel K can be put into 
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the boundary condition. We can treat the + and — quantities as if they were independent. Then, we can derive the 
Feynman rules according to Sctp- The main difference between the tree- level vertices obtained from Sqtp and from 
S lies in that there are negative-type vertices (all time arguments are on the — branch) besides the ordinary positive 
ones. The vertices of the negative and positive types are the same, except that they have opposite sign. At the tree 
level, there is also no vertex of mixed type with both positive and negative time arguments. 

The 2-point GF for Q is defined by G — (TpQ(xi)Q(x2)}, where Tp indicates a path-ordered product along the 
CTP. For simplicity we suppress the Lorentz and colour indices of Q and restore them when necessary. There are 
four types of 2-point GFs characterized by the time-branch (+ or — ). This can be explicitly written as a matrix: 



G = 



(TQ(xi)Q(x 2 )) (Q(x 2 )Q(x 1 )} 
{Q{x 1 )Q{x 2 )) (T*Q{ Xl )Q{x 2 )) 

G++ G+-\_f G F G< 
G-+ G— ) \G> G F 

where T* denotes the counter-time-ordered product; G ++ (G F ) is the ordinary GF and its time arguments t\ and t 2 
are on the positive time branch; G~~ {G F ) is the counter-time-ordered or anticausal GF, with both time arguments 
on the negative time branch; G" 1 (G < ) and G h (G > ) are correlation functions with time arguments on different 
branches. The SE has a similar form: 



, u F n< , 

n=( n>r H- (3D 



The GF and the SE can be expressed in the so-called physical representation by using the following unitary transfor- 
mation: 



U 





n< 


n> 


u F 


G F 


G< 


G> 


G F 



n A o 



where 



^viU i )> u =vil-i i 1 • (33) 



Writing the transformation (B2ft explicitly, we obtain the following relations: 



G A 


= G F 


-G> 


= G< 


-G F , 


IL A 


= n F - 


hn> = 


-n< 


-n T 


G R 


= G F 


-G< 


= G> 


-G T , 




= n F - 


hn< = 


-n> 


-n T 


G c 


= G F 


+ G F 


= G< 


+ G>, 


n c 


= n F - 


fn F = 


-n< 


-n> 



(34) 

with A, R and G denoting the advanced, retarded and homogeneous GF or SF. Note that there is an additional 
negative sign for n >,<; with respect to G > < . The reason is as follows: the SE tensor II is a 1-PI 2-point GF with 
external legs amputated, and the two time arguments of II >,< are on different CTP branches. We also know that 
negative and positive type vertices differ in sign. Therefore, II >:< has an additional negative sign relative to G >,< . 
Such a case, however, does not occur for II F , because both of its time arguments are on the negative branch; they 
thus contribute with the same negative signs, which cancel. 

Keeping these Feynman rules in mind, we can write the DSE (pp[) and (pT|) in the ordinary representation as 



V + II+ \ (G F G<\ 

o -2?-n+ J h {xi) {G> g f J {xi ' X2) 



igfj.xS a c5(xi - X 2 ) 



1 

I 



[ tt> ttF~ I (xi,x)[ r ^ -p {x ,x 2 ) , (35) 



n> n F J ' u >\ g> g f , , 
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G F G<\ ( x>t + n+ \ . . 

G> G F ( x ^){ o -pt_ n + J Xh (X2) 



tl , f G F G<\ ,.fu F n< 



d ( G > g f j (^l.O ( n> h f ) (a') 2:2) , (36) 



or, in the physical representation, as 



V + U+ \ , . / G R G c \ , , 
-V-IL+ ) ( G A ) {X1 ' X2) 



1 n n J C^i,^') (A nc ) (x',x 2 ), (37) 



J G R G c , x . 



g a \ , /©t + n+ \ 

G c G * ) (*i,* a )( Q _ pt _n+ J Xh (X2) 

1 



«3/ J A'5ac^(a;i - x 2 ) 



1 



^'(g-r g c ) fci'*') (n- 4 ^ ) ( x ', x 2)- ( 38 ) 

In Eqs. ( |35f]38| ) , (xi) and (11,2:2) stand for space-time arguments for the functions in front of them. The Lorentz and 
colour indices are written as subscripts. As an alternative way of presenting the above equations, the local term ±II + 
can be absorbed into the definition of H F / F . In the Feynman gauge (a = 1) the differential operators V and Z>t are 
defined as: 



Z>£ = g P „D^[A]D a j[A}-D h a a [A}Df[A) 
+ -D h p a [A]D a a l [A]+gf hm F^[A] 
= g p(T D h "[A}D°?[A} + 2gf hai Fp a [A] , (39) 



and 

= 9 P aD^ ha [A}D^[A}+2gf hm F; a [A}, (40) 

where £)^ a [A] = d IJi S ha + gf' A „ is the covariant derivative in the adjoint representation and Djj ha =d^ 8 ha — 
gjhba^b jg fti e conjugate covariant derivative, where the differential operator acts to the left. We note that Eqs. 
(HHH) are independent of the gauge parameter a, owing to the gauge conditions: D' 1 ^ [A(xi)]Gj^(xi, x 2 ) = and 
G%(x 1 ,x 2 )D^ k [A(x 2 )} = 0. 

Some comments about Eqs. ( p5| - |40| ) are in order. First, we recall that Eqs. ( p9[|40| ) come from the second line of 
Eq. (jl9|). We write it down explicitly: 

*S)iJ>>f) = r^(Q 2 ) + rW mp (AQ 2 )A m , + ir^ n , mp (A 2 Q 2 )A m ,A„, 

= {g pa D h «[A{x)]D a ;[A{x)] + 2gf ha *F; a {A(x)}}s i (x - y) , (41) 

where labels m and p stand for a group of indices: m = (x, p, h) and p = (y, a, i). Multiplying Eq. ( fH] ) by G^(y, z) 
and integrating over y, Eq. (^) gives 

d VG^;^(x,y)G^(y,z) 
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{g p(T D^[A(x)}D^[A(x)]+2gf hai F^[A(x)}} j d 4 y6 4 (x - y)G%(y, z) 
{g p(T D* a [A(x)}D™[A(x)} + 2gf hai F^[A(x)}}G i J i (x ) z) 



:V%(x)G%(x,z) 



(42) 



Multiplying Eq. ( p~l| ) by G| p (z,x) and integrating over x, it becomes: 



d 4 xGf(z,x)iG (0 \: hl a {x,y) 



d 4 xGl h p (z,x){g pa D'^[A(x)}D a ;[A(x)} + 2gf hal F^[A(x)]}d 4 (x - y) 
= j d 4 xG^(z,x){g pa D^[A(x)}D^A^ 

= Gf p (z,y){g p ^ ha [A(y)]D^[A(y)} + 2gf hai F^[A(y)]} 
= G\ h p {z,y)V^{y), 



(43) 



where the differential operator d acting on the <5-function in the second line is changed to — d in the third line. We 



notice that G,}J % (x,y) in Eq. (Ell) should be understood as G ( S' ni (x + ,y + ), where x + and y + are on the positive 



l;hi 



(0);p<t 



(0);p^ 



time branch. In addition, since G ^'.^{x , y) is associated with vertices T^(Q 2 ), r^ mp (AQ 2 ) and ^} n , mp (A 2 Q 2 ), 



Gj j!^(x_, y_) should have an additional negative sign. This is the reason why there is a negative sign before T> and 
Dt in Eqs. (HH). 

Of particular relevance to our further discussion are the specific matrix elements originating from Eqs. ( |37| ) and 
(|38|). The equations corresponding to the upper-right element of (|37| ) and to the lower-left element of (|3q ) are: 



V{ Xl ) +IL + {x 1 ) G c (x u x 2 ) 



d 4 x' 



U c ( Xl ,x')G A (x', x 2 ) + U R ( Xl ,x')G c (x', x 2 ) 



G c (x 1 ,x 2 ) V\x 2 ) + n+(x 2 ) 



d 4 x' 



G R (x 1 ,x')U c (x', x 2 ) + G c (x 1 ,x')n A (x / , x 2 ) 



(44) 



(45) 



where, for simplicity, we suppress colour and Lorentz indices in G and II, and in particular we regard them as matrices 
in colour space. We note that the local SE tensor II + can be absorbed into II fl and H A to make Eqs. (44) and ( |45| ) 
more compact. 



VI. DISCUSSIONS OF THE NON-LOCAL SOURCE 

So far we have not discussed the role of the non-local source term K(A±,Q±) in the DSE. The reason is that 
K(A±,Q±) is non-zero only at the initial time and hence can be put to the boundary condition. To illustrate the 
role of K , we consider a simple example of a non-interacting massive scalar field. The classical Lagrangian density is 
£ = \d p <f)d p (j) — im 2 2 . In the CTP approach the role of the initial density matrix is taken by the non-local source 
kernel K(<f>±). K((f>±) can be expanded into a functional Taylor series in c/>. The linear term K a (f> a of this series can be 
absorbed into the source term J a 4> a - The lowest order term, which does have an effect on the dynamics, is the square 
term -R'ah'^a^b- In this heuristic argument, we neglect higher order terms in this expansion. Using this approximation 
the DSE has a very simple form: 



S 4 (x - y){-d 2 v - to 2 ) + K(x, y) r,A(y, z) = S 4 (x - z) , 



(46) 



where 77 = diag(l, —1); «A(y, z) is a 2-point CTP-form GF which is a 2 x 2 matrix and K (x, y) is a CTP-form matrix. 
The solution of the above equation can be written as 

A(x,y) = A^(x,y) - A^(x,u 1 )r,K(u 1 ,u 2 )v^ Hu2,y) 
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+A<°> (x, u 1 )t 1 K(u 1 ,u 2 )v^ {0} (tta, "3)r?iC(M3, »4)ijA (0) (u 4 , y) 



where A(°) (x, y) is a solution of the inhomogeneous equation 



(-d 2 -m 2 )r]A^(x,y) = S 4 (x-y) 



(47) 



(48) 



which we call a full propagator. The propagator A^°'(x,y) can be written as the sum of the homogeneous solution 

and an inhomogeneous one: A^ ^ = A^„ + A^ m , where *A^P(x ~ v) 1S * ne Feynman propagator and A^ m (x — y) is 
the solution of the homogeneous equation 



{-d 2 - m 2 )r]A hom (x, y) = . 
In momentum space, these propagators become 

A(°>(p) 

where A^ (p) and A^ m (p) are given by 



A ( °\p) + A^ m (p), 



(49) 



(50) 



A (0) (v) 

^hom \P) 



p 2 —ni 2j rie 







-1 



-2m5{p 2 - m 2 ) 

\ U \P 

-2^(5(p 2 -m 2 )[Ai+ ff (0) (p)A 2 



.g(°)(p) 9(- p °)+g(°)( p ) 
0(P O )+9 {O) (P) 9 {0) (P) 



with 



Ai = 



9{-p°) 
9{p°) 



A 2 = 



1 1 
1 1 



In our further discussion we assume that the kernel K(x x ,x 2 ) is given in the following form: 

1 



(2tt)-- 



-^( Xl -x 2 )«5(x?-4)5(x°-g, 



(51) 



(52) 



(53) 



where i,j = ± and anc * are ' respectively, the starting and ending points of the CTP. The kernel K %3 (xi, x 2 ) is 
translationally invariant in space, and its Fourier transform reads 



i^(ki} 



d 4 x 1 d 4 x 2 e lklXl e- lk2X2 K lj (x 1 ,x 2 ) 

^'(koa^ki - k 2 ) , 



(54) 



where we use the same symbol if to denote the kernel in the coordinate and in the momentum space. 

Let us now consider the full propagators A^(x, u x ) and A(°)(u„,x) appearing in Eq. ( |4^ ) at the leading and at 
the end of each term. The values of u\ and are or t^. If one assumes that ijj" — > — oo + and ig = — oo _ then 
A(°);++(x,u 1 ) = A(°) ;> (x,wi) since x° > t+ (x° is definite). Similarly, we have: A(°) ;+_ (x, u x ) = A^'<(x, u x ), 

A(°);-+(x,Ui) = A(°) ;> (x,wi) and A(°) ; — (x,u x ) = A(°) ;< (x, u x ). The ++, H — , — h and components of 

A(°)(it n ,x) are of type <, <, > and >, respectively. In momentum space we have 



A^(p) = -2^% 2 -m 2 )[A 3 +. 9 (0 '(p)A 2 
A<°>( ? ) = -2wiS(q 2 -m 2 )\\ A + g^(q)X 2 



(55) 
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where the subscripts h and e denote the leading and the end propagator, respectively, and A3 and A4 are defined as 



{ e(p°) o(-p°) ) 



A 4 = 



o(-q°) e(- q °) 

0(q°) 6(q°) 



(56) 



The propagators A°(uj, Ui+i) appearing in Eq. (^) have both of their time variables pinched at to. Thus, they are not 
bound by the above arguments and their four components are equal. Finally, in momentum space, the full propagator 
(Et]) becomes 



A(p, q) = A^(p)S 4 (p -q)- A( 0) (p)/C(p)A(°)( 9 )^(p - q) 

where /C(p) is given by 

AC(p) = V K(p) V ~ V K(p)7 1 A^(p)\ 2 r 1 K(p) V 

+f 1 K(p)f 1 A (0) (p)X2T]K(p)f ] A (0) (p)X2T]K(p)f ] - ■ ■ 

and A(°)(p) = / dp°A^(p). Collecting the sum in the above equation, we obtain: 

/C++(p) = c[(l + AB 2 )K++(p)+AB 1 K+-{p) 

fC+-(p) = c[-{l + AB 1 )K+-{p)-AB 2 K + +{p) 

K,-+(p) = c\-(l + AB 2 )K- + {p)-AB 1 K—(p) 

K—(p) = c\(l + AB 1 )K—(p)+AB 2 K- + (p) 



(57) 



(58) 



(59) 



where A = A< )(p), B x = K++{p) - K~ + {p), B 2 = K (p) - K+-(p) and C = 1/[1 + A(B 1 + B 2 )}. 

There are three types of contributions to the second term of Eq. (p7|): A3/CA4, A3/CA2 + A2/CA4, and \ 2 K.\ 2 . We 
denote them as ii, I 2 and ^3 respectively: 



where Iq is defined by 



h = I [IC ++ S(p +q )+IC + -5(p - q°)}\ 2 , 

h = 7 .9 (0) (P)(^ ++ + >C + -)[6(P - q") + S(p° + q°)]\ 2 , 
h = W 0) (P)] 2 (/C++ + JC + -)[S{p° - q°) + S(p° + q°)]X 2 , 



h = Sir 2 5{p 2 -mf)^ (p-q) 



(60) 



2E r , 



(61) 



with E p = ^/p 2 



In Eq. (tefj) we have used: K ++ = K and K + = K + , which implies that JC ++ = JC 



and JC + = JC + . We have also used the following formula 



S(p 2 - m 2 )S(q 2 - m 2 )5 6 {p - q) = 6{p 2 - m 2 ) — [5(p {> - q") + S(p a + q ^ 

2E p 

With Eq. @, Eq. @ becomes 

A(p,q) = A ( -°\p)5\p-q) + h+I 2 + h 
M 



(62) 



= AZ'(p)-2mS(p 2 ~m 2 ) A 1+ /(p)A 2 \5\p-q) 



-2mS(p 2 - m 2 )/'(p)A 2 5(p + 9o)<* 3 (p - q) , 



(63) 



where f(p) and /'(p) are given by 



2tt 



/(P) = S (0) (p)+*tH /C+-(p)+ ff (0) (p) + ( 5 (0) (p)) 2 /C++(p)+/C+-(p) 



/'(P) 



. 2vr 



/C++(p)+ .g (0) (p) + (.9 (0) (p)) 2 /C++(p)+/C+-(p) 



-,(°) 



(64) 
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The appearance of 5(p° + q°) in Eq. ( |63| ) could, in general, result in the time dependence of the full propagator. 
However, if we assume that the initial time to is in the remote past (i = — oo), while x° and y° of the full propagator 
are finite, we can consequently drop 5(p° + q°). This is because, when going to coordinate space, this term would 
generate a factor e~ lp ( x +v ~ 2t ' that vanishes according to the Riemann theorem. Under this assumption the full 
propagator in Eq. ( p3[ ) can be finally written as 

A(p, q) = | A#> (p) - 2m5(p 2 -m 2 )[\ 1+ /(p) A 2 ] | 5 4 (p - q) , (65) 

which has the standard form generally assumed in the literature. 

In this section we have discussed the effect of the non-local source kernel K on the solution of the DSE in a 
simple free scalar field model. In terms of this model we have derived the full propagator and shown that K entering 
the homogeneous part of the solution brings its time dependence and breaks the time translational symmetry. If one 
assumes, however, that the initial time is in the remote past, one then finds that the time dependence can be neglected 
and thus the time translational symmetry is restored. In this case the effect of K can be collected to the homogeneous 
solution in such a way that K only corrects the distribution function. In equilibrium, this distribution function is just 
the Bose-Einstein distribution. In the above free scalar field model that neglects interaction the equilibrium cannot 
be, however, reached. For a more complicated cases such as QCD, the effect of K on the solution is more involved. 
Generally K can be treated as a special kind of boundary conditions pl| . Therefore, the kernel K is absent in the 
derivation of the transport equation that is presented in the following sections. 

The non-local source K is defined as the matrix element of the density matrix p on the initial states 4>i,2 

(0i,Ok# 2 ,O) = exp{iK(<p±)] . (66) 
The kernel K can be expanded functionally as follows: 

K{(j)±) = K + J d 4 xK a (x)<f> a (x) + ^ J d A xd 4 x'K ab a {x)(f> b (x') + --- (67) 

with a, b = +, — . In general K is a complex functional of the fields. From Eq. (^3|) and (|65|) the coefficient K, of the 
quadratic term in the functional expansion of the source K should be purely imaginary to ensure that the spectral 
function is real. For a general form of the density matrix p ~ exp(— J d 3 k(3k(^ k a\ t ), one can indeed check that when 
expanding the kernel K up to the quadratic term the coefficients are all imaginary [see Eq. (2.31) 0]]. 

It is interesting to note that the situation here is similar to the pinch singularity, which arises when the time 
variation of the distribution function |35| is neglected. The possible connection between the non-local source kernel 
and the pinch singularity will be discussed elsewhere. 



VII. KINETIC PART OF THE TRANSPORT EQUATION 



In Section V we derived the DSE for a gluon plasma in the CTP formalism. The resulting DSE summarized in Eqs. 
(||) and (f|| is a non-linear integro-differential equation, which cannot be solved without further approximations. 
The essential approximation usually made in the literature pq| is based on the two-scale nature of high energy QCD. 
There are two typical scales for a multiparton system: the quantum scale, which characterizes quantum fluctuations 
or parton self-interactions, and the statistical-kinetic scale, which measures the range of interactions between quasi- 
particles. These interactions may be described in a semiclassical way if these two scales are well separated, i.e. if the 
local density of quasi-particles is smaller than a critical density where particles begin to overlap. The above situation 
is well suited to the case of ultra-relativistic heavy-ion collisions. Shortly after two highly Lorentz-contracted nuclei 
pass through each other, a very strong background field is formed, followed by the production of very high energy 
partons. Since this occurs very early and at a very short space-time scale, it is purely a quantum process, thus 
the quasi-particle-based semiclassical or the kinetic description generally fail. As time goes on, the dense system of 
partons undergoes an expansion and the local density may fall down to a level where the quantum and classical scales 
can be well separated. Through multiple collisions, the parton system can thermalize and then its bulk properties 
can be described in terms of hydrodynamics. 

With the above physical scenario in mind, let us assume that the space-time is discretized into cells of a size chosen 
such that the separation between the quantum and kinetic scales [fl3| [Ji] is optimized. Then, the correlation between 
different cells will be negligible. The 2-point correlation will not vanish, only when two space-time points lie in the 
same cell. Consequently, in the multiparton system and in a given cell, one can neglect spatial inhomogeneity of the 
local gluon and quark densities. Within each cell, one may therefore describe the short-distance quantum dynamics 
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analogously, as in vacuum or in a homogeneous medium. The inhomogeneity of the spatial parton distribution 
associated with particle collisions appears only when moving from cell to cell. 

Let us define a mass scale /i as the separating point of the quantum and the kinetic scale. This implies that one may 
characterize the dynamical evolution of the parton system by a short-range quantum scale X qua < ~ , and a long-range 

kinetic scale Xkin > ^ ■ The low-momentum collective excitations that may develop at the particular momentum scale 
gH are thus well separated from the typical hard gluon momenta k > fx, provided that g <C 1- The effect of the 
classical field A on the hard quanta involves the coupling gA to the hard propagator; it is thus of the order of the soft 
wavelength ~ l/(gfi). Hence, we have the following characteristic scales: 



y = xi - x 2 



1 

A* 



1 



d v = — d 2 ) ~ fi 



X = -(xi + x 2 ) 



1 



dx = di + d 2 



9V 



gA(X)~g^, gF[A(X)\ 



2 2 

9 V 



(68) 



where X is the central point and y the difference of two coordinates x\ and x 2 in a 2-point GF. We see that X labels 
the macroscopic kinetic motion whereas y characterizes the microscopic quantum distance. 

In order to take advantage of the assumed separation of scales we first express all 2-point GFs appearing in Eqs. 
(|35|), ( |36| ) and (|37j), ( |38| ) in terms of the new variables X and y. Then, to derive the transport equation from the DSE, 
we perform a gradient expansion of these GFs under the conditions (|68|). From Eqs. ( f44| ) and (45|) it is clear that we 
will deal with such an integral as I = J d A x'Ii{xi : x')G(x' , x 2 ). In terms of X, y and y' = x' — x 2 coordinates, the 
integral and its Fourier transforms, with respect to the relative distance y, can be expressed in the gradient expansion 



/ = / d 4 y' U(X, y - y')G(X, y') + -8 X U(X, y - y')y'G(X, y') 



F[I] 



-(y-y')H(X,y-y')dxG(X,y') 
d 4 y e igy j = n(x, q)G{X, q)+ l - 



d q U(X, q) ■ d x G{X, q) 



~d x Tl(X,q)-d q G(X,q) . 
Similarly, the Fourier transforms for I\ = H{x\)G{xi 1 x 2 ) and I 2 = G{x\ 1 x 2 )H{x 2 ) are 

F[h] = U(X)G(X,q) - l -d x ll(X)-d q G(X,q) , 



(69) 



F[I 2 ] = G(X,q)U(X) + -d q G{X,q) ■ dxll(X) 



(70) 

The 2-point GFs in Eqs. (|35j), ( |36| ) and (^7(), (|38| ) are not gauge covariant. In order to obtain a gauge-covariant 
transport equation, we must use a gauge-covariant 2-point GF defined by 



G{X,y) = V{X,x 1 )G(x 1 ,x 2 )V(x 2 ,X) , 
where V(zi, z 2 ) is a Wilson link with respect to the classical background field given by 



V(zi,Z2) = T P exp ^ig J dz^A^ 



(71) 



(72) 



where the integral stands for a path integral from point z 2 to z\ and Tp denotes the ordered product along the path; 
note that the path here is defined in coordinate space. The Wilson link V(zi, z 2 ) transforms as 

V(zi,z 2 ) -> U( Zl )V( Zl , z 2 )XJ-\z 2 ) , (73) 

where U(z) = exp(igu> a (z)t a A ) is the gauge transformation under which the GF G{x\,x 2 ) transforms as 

G( Zl ,z 2 ) - U(z 1 )G(z 1 , z 2 )U T {z 2 ) = U( Zl )G( Zl , z 2 )U-\z 2 ) , (74) 

which involves transformations at two different space-time points. However, the CGF G(X,y) transforms as 

G(X, y) -> U(X)G(X, y)U-\X) , (75) 
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where only the transformation at a single point X is relevant. The gauge-covariant Wigner function is the Fourier 
transform of G{X, y) with respect to y: 

G(X,q) = J d 4 yG(X, y y™ . (76) 

Obviously G(X,q) transforms in the same way as G(X,y) according to Eq. (|75|). 

In general, the integration path in the Wilson link may be of arbitrary shape, provided that its two end points are 
fixed; the gauge-covariant Wigner function is thus not uniquely defined. This ambiguity can be removed by requiring 
that the Fourier spectrum variable q in the Wigner function G(X, q) corresponds to the kinetic momentum in the 
classical limit. This constraint is fulfilled by the straight-line path p9j. Therefore, in all our future calculations we 
imply the straight-line path. The link operator defined on the straight-line path has the following properties: 

V(z 1 ,z 2 )V(z 2 ,z 3 ) = V(zx,z 3 ) , 
V(z 1 ,z 2 )V(z 2 ,z 1 ) = 1 , 

V f ( Zl ,z 2 ) = V- 1 (z 1 ,z 2 ) = V(z2,zi) . (77) 

Obviously, the Wilson link is unitary since A = A* and forms a group. With the above properties, we can write the 
inverse relation for Eq. ((n|) as 

G( Xl ,X2) = V(x u X)G(X,y)V(X,x 2 ) . (78) 

Consider a straight-line path from z 2 to z\ described by the equation z(s) = z 2 + [z\ — z 2 )s with s = [0, 1]. The 
variation of the path characterized by small changes of their end points dz\ and dz 2 causes the following variation of 
V" : 

5V(zi,z 2 ) = igA{zi)dziV{zi,z 2 ) - igV{zi,z 2 )A{z 2 )dz 2 
-ig dsv[z u z{s)^F^[z{s)^v[z{s),z 2 ^ 

■(zi - z 2 )^ dz 2 + (dzi - dz 2 )s . (79) 



Using the above equation, we obtain 



d xlfl V(x u X) = --igV{x 1 ,X)A^{X)+igA fl {x 1 )V(x 1 ,X) 



1 

— 2 

■ 3 r. / 3 1 

-tg-VuFu^-xi + -X 2 

d Xl ^V(X,x 2 ) = ^igA^X)V(X,x 2 )-ig^y v F u ^ Xl + ^x 2 ) . (80) 

The first equation can be written also as 

D x ^V(x 1 ,X) = ~igV(x 1 ,X)A^X)-ig~y v F„^x 1 + ~x 2 ) , (81) 

where the l.h.s. is O(gn) and the two terms on the r.h.s. are O(gfi) and 0(g 2 /j,) respectively. 
Taking the Hermitian conjugate of Eq. (|8fj| ) and interchanging x\ and x 2 , we obtain: 

d X2li V(X,x 2 ) = UgA^X)V{X,x 2 )-igV{X,x 2 )A^x 2 ) 

■ 3 r, ( l 3 \ 

-W-^VuFu^y-X! + -x 2 J , 

d xni V(x 1 ,X) = -ligV(x 1 ,X)A li (X)-ig^y l/ F 1/fi (^x 1 + ^x 2 y (82) 



From the above results we can also find that 



D XlV {VxGV 2 ) = ^(DxvG^ + V^dyvGW-ig^yxFxvGVi 
1 „ x . „ 1 



2 igViGA u V 2 - ig—ViGy\F\ v , (83) 
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with the following compact notations: Vi = V(xx,X), V 2 = V(X, x 2 ), G = G(X,y) and D X1U = D l/ [A(x%)]. The 
terms on the r.h.s. of Eq. (83) are of orders g/x, /i, g 2 fi, g\i and g 2 /i, respectively. We also make the approximation 
F\ v (\x\ + §£2) ~ F\ V (X) in all terms involving F\ v in the above equation, because the corrections are of 0(g 3 fj,). 



From Eq. fl83| ) and its Hermitian conjugate (ViGV 2 )D\. 2V we can derive the gauge condition for the gauge-covariant 
GF. The background gauge condition is given by D® b [A(x)]Q b ^(x) = 0, thus the gauge condition for the Green function 
with respect to X\ reads 



D?[A(x l )]G%,(x l ,x 2 ) = 0, 
where G is one of G > , G < or G c . With respect to X2, the gauge condition is 



D^[A(x 2 )}G b l(x u x 2 ) = 



whose complex conjugate is given by 



6 ac d 2 „ +ig[(t d A rrK(x 2 ) G^(xi,x 2 ) = 



(84) 
(85) 
(86) 



where we have used the fact that A and G are real. The SU (3) generators in the adjoint representation are Hermitian 
(m.) = t% therefore we have 



5 ac d 2v +tg(t d A ) ca Ai(x 2 )\G%(x 1 ,x 2 ) 
--G%(x 1 ,x 2 )D$> ca [A(x 2 )}=0. 



(87) 



From Eqs.@ 



, and (B7p we find that 
D^iVxG^Vv) 



f 1 - 3 ~ 

VA -{Dx^G^u) + {dynGpv) - ig-y\F x ^G^ v 

+~igG IMU A tl - ig^G^y x F x Av 2 = , 



and 



(ViG Vfl V 2 )Dt 2fl = v^(G^D x J-(d y »G Vfl )-ig~G^y x Fx 
-^ig^G vt , - ig^y\F XtJ .G u Av 2 = , 



(89) 



where we kept only terms up to 0(g 2 fi). Taking the sum and the difference of the above two equations, we derive the 
gauge condition for the gauge-covariant G: 



-^dx^G^y + dy^G^] + -ig G 
-zdx^G^ +d yll G{ llv -) + -ig G^^A^ 



0, 

= 0, 



(90) 



where Gs^, v \ and Gy^ are defined by G^ v + G ufJi and G^ v — G vil respectively. For simplicity we also neglect the terms 
of order higher than 0(g/i). 

If we assume that G is symmetric in its Lorentz indices, then up to 0{gn) we obtain 



dx^G^v + ig[Gfj, v , A^] — , d yil G^ v — . 

The second equation tells us that G^ v is transversal up to O(gfi). 

Taking the second covariant derivative of Eq. (|83|), we obtain the covariant d'Alembertian operator: 



(91) 



DIJVL&V2) = -ViiDxtyVz + Vridy-DxG^ + ViidlG^ 

3 3 
-ig^y\F\v{D X vG)V 2 - ig-y\F X v{d yu G)V 2 

+ig\vi{D x „G)A u V 2 + ig^ViG(dxMV 2 

-ig^VxiDxuG^xFxv - ig^Vi{dy V G)yxF >v 

+igV 1 (d yv G)A l/ V 2 - X -g 2 V x GA 2 V 2 , 



(92) 
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where we kept only terms up to 0(g 3 fi 2 ), neglecting all higher order ones. The resulting equation for (ViGV 2 )D^ 
can be derived by taking the Hermitian conjugate of D^. (V± GV 2 ) and then interchanging x\ with x 2 . 

{V x GV2)Dll = ^V 1 (GD x i )V 2 -V 1 (Gd y -D x )V 2 + V 1 (d 2 G)V 2 
-i9^V 1 {GD^ Xu )y x F Xv + ig^dy^yxFx, 
-igh/ x A v {GD Xv )V 2 - ig~Vx{d Xv A v )GV 2 
-ig\yxFx,{GD Xu )V 2 + ig\yxF Xv {d y „G)V 2 

+igV 1 A u (d yv G)V 2 - - A g 2 V 1 A 2 GV 2 . (93) 



We take the difference between Eq. ( |92D and Eq. ( |93|) : 

V 1 AV 2 = D 2 X1 (ViGVi) - {VxGV 2 )DV , (94) 

where A is defined by 

A = d y ■ D X G + Gd v -D x - ig^yxF Xu (D Xv G) + ig^(GD Xu )yxF Xv 
-igyxFxv{d V vG) - ig{d yv G)y x Fx v 
-ig~(D X uG)yxFxv + ig^y x F Xv (GD Xv ) 

+ig(d yu G)A v - igA v {dy„G) . (95) 
Keeping terms up to 0(g 2 ^, 2 ), we obtain 

A = d y ■ D X G + Gd y -D x - igyxFxv(dyvG) 

-ig(d y „G)y x Fxv + ig{d yv G)A v - igA v {d yv G) . (96) 

Taking the difference of the l.h.s. of Eqs. (Q and (^|), and using Eqs. (|j|) an d (H) we get 

A' = dy D x G c ai + G% d y -D x + igyxFvxidyvG^) + ig(d yu G^)y x F v x 

+ig(d yu G c ai )A v - igA v {d yv G c ai ) ~ 2ig{F a pG% 1 - G c ap F Pl ) , (97) 

where we restored the Lorentz index for G (X, y) and suppressed two Wilson links in front of and behind the above 
expression. With respect to the relative coordinate y, the Fourier transform reads 

F[A'} = -i[q-D x G% + G c ai q-D x +g qu F v x{d q xG%)+g{d q xG%)q v F v x 

+ig{G c ^q -A-q- AG%) + 2g{F aP G% - G%F Pl )) , (98) 

where G c = G C (X, q), A = A(X) and F = F(X). 

Equations (J^7§ and ( |98|) just describe the kinetic part of the transport equation. The collision part will be derived 
in the next section. Neglecting the collision terms, which are of order at least g ' fi , we have the kinetic equation in 
the following compact form: 



q ■ d x G^ + ig(G^q -A-q- AG^) 
+ ^9QuF uX {d q xG^) + ^g(dqxG% y )q v F v \ 
-9{F aP G c M - G%F Pl ) = . (99) 



The above equation is located at the collective coordinate X and is gauge covariant under the local gauge transfor- 
mation U(X), i.e. it transforms as U{- ■ -)?7 -1 . Indeed, noting that both F^ v and G^ 7 are gauge covariant and d q ^ 
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does not affect U, it is obvious that the last three terms are gauge covariant. To verify that the first two terms are 
also preserving gauge covariance, we explicitly write down their transformations: 

q-d x G% -» Uiq-dxG^U- 1 + (q-d x U)Gg ( U- 1 

+UG%q ■ OxU- 1 , 
igG^q ■ A -> igUG c ai q ■ All' 1 - UG%q ■ dxU' 1 , 
-igq-AG^ -► —igllq- AG^U^ 1 

-(q-dxU^U- 1 , 

from which it is clearly seen that the sum of the first two terms in Eq. ( p9| ) indeed transforms as U(- ■ -)C/ _1 and 
therefore preserves the gauge covariance. 

The quantum kinetic equation ( |99"| ) was derived in a quite general and transparent manner in the context of BG- 
QCD and CTP formalism. No further approximations or requirements going beyond gradient expansion were used 
to obtain Eq. (9^). We note that a result similar to Eq. ( p9] ) was previously obtained in Ref. j3(|. There, however, 
based on Ref. the transport equation was derived by making the gradient expansion of the equation of motion 
for the Wigner function (not in CTP formalism). Additionally the author of Ref. |36| made the derivation in the 
fundamental color space (not the adjoint space) in QCD (not in BG-QCD). Finally, in Ref. |56), Eq. j99| ) was obtained 
by assuming that the Wigner function is proportional to the quadratic product of the generators of the fundamental 
representation. Our approach is quite general and does not require any specific assumptions on the structure of the 
Wigner function. 



In the following we show that Eq. (99) is a natural quantum generalization of the classical Boltzmann equation. 
In particular the colour charge precession will be explicitly identified in the quantum description of the colour charge 
kinetics given by Eq. (|9E 



The classical kinetic equation for the colour singlet distribution function f(x,p, Q) is |ll], [ll], 23 [M. 25, p6| : 

p^-gQ a F^d pv 

- 9 r bc Al(x)Q c d Q .}f(x, P ,Q) = , (100) 

where Q a is the cl assic al colour charge and a = 1, — 1. 

Comparing Eq. ( |l00| ) with the quantum expression (Q) , it is clear that the colour singlet distribution function / is 
replaced by the gauge covariant Wigner function G c , which is a colour matrix in the adjoint representation. One can 
also recognize that the first, third and fourth terms of Eq. ( p9[ ) arc the quantum generalization of the first two terms 
in Eq. ( |l 00| ) - The last term in Eq. ( p9[ ) appears from the covariant operators and hence it does not appear in the 
classical equation. This term can be written in different form by using generators of Lorentz transformation in vector 
representation. A similar term can be found for the quark, but expressed through generators in spinor representation. 

Particularly interesting is the appearance of the second term in Eq. (|9^) . We have seen that its presence is crucial to 
assure the gauge covariance of the Vlasov equation. This term has an interesting physical meaning. It is the quantum 
analogue to the color charge precession in the classical kinetic equation. To see this more clearly, one can expand 
G^g(X, q) with respect to the expansion parameter gT a A^(X) from the Wilson link in Eq. (|7l|). This expansion can 
be also understood as the result of the AQQ, AQQQ and AAQQ vertices. Then we have 

G%{X,q) = N QaP (X,q)+gT a A«(X)N la ^(X,q) 
+g 2 T a T b A«(X)AUX)N 2ap .^(X, q) 

+ •••, (101) 

where T a are quantum analogues to the classical color charges Q a ; Ni a p(X,q) with i = 0, 1,2, • • • are color singlet 
functions. Each term of the expansion corresponds to an order of the color inhomogeneity in the gluonic medium due 
to its interaction with the background field. If the background field is refered to the soft mean field, its mag nitud e 



should vanish when the system approaches equilibrium. Then only the first singlet term survives in Eq. ( 101 ), 
which means the color homogeneity of the gluonic medium. This is somewhat similar to the multipole expansion 
for an electromagnetic source where the moments of dipole, quadrupole etc. describe increasing orders of spatial 
inhomogeneity for the elec trom agnetic charges. In weak coupling, as the lowest order approximation, we keep only 
the first two terms in Eq. ( |l0l| ) . Then the second term of Eq. (|9^) becomes 

ig{G%q -A-q- AG%) ~ - 5 r bc ^^T c 9 T <, G% (102) 



which reproduces the classical color precession term in Eq. (10C) 
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The covariant derivative Dx ~ gA(X) ~ g/i and therefore first two terms of Eq. (^9|) are at leading order, 0(g/i 2 ), 
while other terms are at subleading order, 0{g 2 p 2 ). In the vicinity of equilibrium the natural scale in the system is the 
temperature T. The mean distance between particles is of the order of ~ 1/T, while X/(gT) characterizes the scale of 
collective excitations ||^] . For small coupling constant g these two scales are well separated. The covariant Wigner 
functions can be expanded around their equilibrium values: G — + SG, where the equilibrium function G^ is 
a colour singlet and the fluctuation SG ~ g 2 G. Typical scales are q ~ T, Dx ~ g 2 T, gF ~ (Dx) 2 ~ g 4 T 2 . Thus at 



survives and the precession term vanishes because of the color-singlet 
fl99) with respect to 5G corresponds to the equation formulated in the 



leading order, only the first term of Eq. ( 991 ) 
nature of & ' . The linearized version of Eq. 
background Coulomb gauge in Ref . [E8| . 

The quantum fluctuations near equilibrium were also considered in Ref. |37| in the context of the classical collisionless 
transport equation. It is quite natural to carry out the same study from our quantum approach. First, BG-QCD deals 
with the classical field and the quantum fluctuation in a systematic way. The quantum field plays the similar role to 
the field fluctuation in Ref. p7| . Second, in the quantum approach, corresponding to the phase -space distribution, 
we deal with the GFs which can be expanded around their equilibrium values following Eq. (101). One also needs to 
complete the equations by including the field equation (Q) where the averaged induced current is related to the 2- and 
3-point GFs which finally depend on G. Thus it can be expanded according to Eq. ( |l0l|) as well. 

The analogy and differences of the quantum and classical Boltzmann equations can also be clearly exhibited when 
formulating the equations for the colour moments. Corresponding to Eq. (|9^) one gets 



q ■ d x h ai + gq v F^ x d q \h a ai 
+9{F a a ph a Pl - h a a0 F^) = 
q-dxK J +gf abc q-A b h c ai 

+gq,F b ^d q ^K 
+g(F b ht 



Q 



where we define: h aj = Tt(G£ 7 ), h a ai = Tr{T a G%) and h% = Tr{T a T b G^). 
The classical equations for colour moments of f(x,p, Q) are given by 

p ■ d x f(x,p) - gp^ v d pi/ f a {x lP ) = 

P -d x r(x,p)+ 9 r bc p-A b (x)r(x,p) 



9P„F b v d pv f ab (x,p) = 



(103) 
(104) 

(105) 
(106) 



where f ab (x,p) 



f(xj),Q), f(x, 



fdQf ^Q) and f a (x lP ) = / dQQ a f(x,p,Q). 



Comparing Eq. ( |l03| ) with (105) and Eq. ( |l04|) with (|10q ) one sees that, apart from terms like (F a php 1 — h a pF^), 
which come from the covariant operators, t he q uantum and classical equations have a similar structure. The identifi- 
cation of the colour precession term in Eq. (|99| ) is straightforward. 



VIII. COLLISION PART 

In this section we will derive the gauge covariant collision part of the transport equation in a pure gluon plasma. 
The collision part is of the higher order in g as it is suppressed at least by g 2 relative to the kinetic part. 

The collision part is derived by taking the difference of the SE terms in the r.h.s. of Eqs. (Eij) and ( fi"5"| ) and then 
performing a gradient expansion. First, these equations should be expressed in terms of the gauge covariant Wigner 
functions G and EL 

The integral of Eq. (p9) can be written as 



xx, X 

y 



y 



,y-y')v[X 



—,y )V(X ^— ,z 2 



d 4 y'v( 
xG(X- 

V( Xl ,X){ J d*y'[n(x+?L,y-y')G(x 
-ig^-A(X)IL(X,y- y ')G(X,y') 



X 



y-y 



y-y 
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-fl(X,y-y')G(X,y')ig?—y-.A(X) 



+U(X, y - y')ig y - ■ A(X)G(X, y')\ W(X, x 2 ) 



+higher order in y or y' 



(107) 



where we have made a gradient expansion of the Wilson link operators and kept only leading terms in y or y' , applying 
Eq. ©: 



V 



(x u X+ V -) » V(x u X) 



1 - i<?| ■ A(X) + 



,-, v y-y 

Vl A — ,x 2 



V 



X- 



y-y 



l+Z ff |-A(X) . 



(108) 



The Fourier transform of the term inside the curly bracket {■ ■ •} in Eq. (107) reads 

F[I] = ti(X, q)G(X, q) + % -[ [d q fl(X, q)] ■ [G(X, q)D x ] 

+[d q U(X,q)] ■ [D x G(X,q)} - [D x U(X,q)} ■ [d q G(X,q)} 
-[n(X,q)D x ] ■ [d q G{X,q)]) + ^gd q [n(X 7 q)A(X)G(X,q)} 



d„U ■ GDI + <9„n • D X G - D X U ■ d„G ~ UDl ■ d Q G 



-gd g (ilAG) 



(109) 



where in the second equality we suppressed the arguments of G(X, q), H(X, q) and A(X). In the transformation from 



the coordinate to the momentum space we use the following replacement: d y 



-iq and y 



-id„ 



In the above equations all terms containing the derivatives d q and dx are suppressed by g relative to the leading 



term IIG. Keeping only the lowest order contribution to Eq. (109), the collision term reads 

Icon - G R U C + G C U A - U C G A - tl R G c , 



(110) 



where we drop the Wilson links V(xi,X) and V(X, x 2 ) as they can be cancelled with those in the kinetic part. To 
further simplify the collisions term, we use the following relations: 



G R = \{G A + G R ) + \{G> - G<) / U R = Ufl A + tl R ) + i(fi< - ft 



G A = \{G A + G R ) — \{G y — G < ) ' \n A = i(n A 



ii R ) - i(n< -n>) 



(in) 



to obtain 



tcoll 



G R fl c + G c fl A - fl c G A - fl R G c 



= {G<,n>}-{G>,n<} + - 



G> +G < ,fl A + U R 



< r<A 



n > + -n>,G 



G 



(112) 



where G and II are the matrices in the colour and Lorentz indices. 

While the transport equation is derived by taking the difference of Eqs. 
gives the mass-shell equation: 



2) and (Eq) , the sum of the two equations 



1 



1 



Q-ldx)G + -ig (8 X -A),G 



-/{a-a,g} 



~g 2 A„GA„-^ig 



(d Xv G),A v 



F\,v, (d q \G) 



ig 



= I (g r u c + G c fl A + fl c G A + fl R G c 
This equation provides the constraint on the gauge covariant GF. 



(113) 
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The collision term (110), (112) combined with the kinetic part ( |99| ) gives a complete result for the transport equation 
for a pure gluon plasma. A possible way of obtaining some physical insight and interpretation of this equation is to 
consider the simple case that the system is near equilibrium. In this case, we decompose the GF G and the SE tensor 
f± as: G = G(°) + SG and ft = 11^ + 511, where their equilibrium values G^ and LT^ are colour singlets and SG and 
511 denote deviations from equilibrium. Obviously, both the collision and kinetic parts vanish for and IT - 1 . 

In the vicinity of equilibrium the natural scale is the temperature T. The mean distance between particles is of the 
order of 1/T, whereas l/(gT) characterizes the scale of the collective excitations [p| p9| . In the weak coupling limit 
for g <C 1 these two scales are well separated. Therefore we have: q - T, D x ~ifT and gF ~ (D x ) 2 ~ g 4 T 2 and 
the fluctuations SG ~ g 2 G^ and 511 ~ g 2 fl^°\ In the leading order the Boltzmann equation reads 



q ■ OxSGaj 



gF XlJ q v d qX G^ c + g(F a pG^ c - G% c F fil ) 



= t 2 



G (0)< Afr> 



OLl a0 U 01 



U a0 OlL 0i 



rfr< M0)> 



+/>r« n (0)> 4-f\ ( -° s>> /ir< 



SG>Il 



(o)< 



a0 s - s -fil 



n 



(0)< 

<X0 



6G 



0i 
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where the color commutator [&°^ c ,q ■ A] has been neglected, since G^ is a colour singlet. In deriving the above 
equation we have also used the following approximations: 



//=V(0)> 



(5G> 



a0 







^)(Gf> + G%<), 



n 



(0)R, 



01 



n 



(0)R 



a.0 



)(SG> 7 



T 0i 

*01> 



SG<) 



(115) 

We recall that up to 0(g 2 T) (here we have fi ~ gT), the gauge condition (|}l|) requires that G > , G < and G c be 



transversal. If we further apply the transversality conditions in the DSE one can verify that Eq. (115) indeed holds. 

Incorporating the transversality condition, we can parametrize the GF, Gf^l G and SG^i C , in the following 
form: 



f a0 (X,q)=T a f 3 f(X,q) 



(116) 



where / stands for G > l < l c , G^ > l < l c or 8G > / < ^ C . We have also assumed that SG C has the same structure as G^ c , 
except that G^ c is a colour sing let. Note that we have separated all Lorentz indices into the transversal projector 
T a = ga0 — Qaqp/q 2 in Eq. ( |116| ). Then we can use Lorentz scalars / to express the Boltzmann equation (note that 
these /'s are different from the one used earlier). 
Under the assumption that 



SG<p(X,q) = SG> (X,q) = -SG%(X,d) , 



(117) 



and inserting Eqs. (116) and (117) into (114), we find 

T 



q ■ d x SG c - gF Xv q v d qX N^ 
g(F a 0Tp~y — T a ffFp 7 )& QS>c 
iT ai {GW<5iL> - &°>6IL<\ 



(0)< 



ill 



(0)> 



(118) 



where we used the following notation: 



2T Q7 <5n> /< = Tapsn^ 
2T Q7 n(° )>/< ee T aP snf>/< 



T Pl » 

xrr(°)>/< rp 
5Il a0 T Pl 



(119) 



One can recognize in Eq. ( |11S| ) that the first term on the r.h.s. is the collision and the second the damping term, the 

damping rate being given by [iflj!^ — ilD ^^] . The physical interpretation and analysis of these terms can be found 
in the recent pap ers by Blaizot and Iancu p8| p9| . 

Equation ( 1 IS ) is a linearized version of the Boltzmann equation in a pure gluon plasma with respect to the off- 
equilibrium fluctuations. The linearized equation was previously derived in Ref. p8[. In our approach, however, this 
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equation was obtained in the covariant background gauge, whereas in Ref. [E8| it was done in the Coulomb background 
gauge. In the Coulomb gauge, the physical polarizations are entirely contained in the spatial gluon propagator and are 
independent of the Coulomb ghost. In the covariant gauge the physical transverse degrees of freedom are mixed in all 
components of the gluon propagator. Hence, the ghost diagrams are necessary to cancel the unphysical polarization 
and to guarantee the unitarity. Much as for the gluon, one also needs to introduce covariant GFs for ghost fields and 
formulate their evolution and transport equations. Note that only the collision part contains contributions from the 
ghost because they appear in the SE diagrams. 

IX. SUMMARY AND CONCLUSIONS 

In this paper we have presented a systematic derivation of the quantum Boltzmann equation for a pure gluon 
plasma. First we have developed a functional method to derive the DSE in the BG-QCD. The 1-PI vertex and the 
CGF were defined by the functional derivatives of their generating functionals with respect to the field average and 
the external source, respectively. The bare vertex was derived by taking the functional derivative of the classical 
action with respect to the corresponding field. 

We have started our derivation by expressing the classical action in the DeWitt notation, which, in our opinion, 
results in a simple structure of the formalism. Then, taking a functional derivative of the action with respect to the 
gluon field, we derived the equation of motion. We recursively used the relations between the generating functional 
for the CGF and that for the 1-PI vertex to express a higher-rank GF in terms of the lower-rank CGFs and 1-PI 
vertices. Using this method, we easily derived the DSE for the 2- and 3-point GFs. The current approach has the 
great advantage that it can treat a non-local and a local source term in the same way and that it can produce all 
needed Feynman diagrams automatically. Hence our method is easy to implement by a computer algorithm that can 
generate all Feynman diagrams for a given process. 

We gave a heuristic discussion of the effects of the non-local source kernel K on the solution of the DSE for a free 
scalar field. The role of K is equivalent to that of the initial density matrix. In the absence of a kernel, the general 
solution of the DSE is the sum of the Feynman propagator, which is the solution of the inhomogeneous DSE, and the 
solution of the homogeneous DSE. The homogeneous solution involves a particle momentum distribution function, 
which is just, in equilibrium, the Bose-Einstein distribution. We showed that, if the initial time is in the remote past, 
the effect of K can be collected into the momentum distribution function. Thus, the structure of the homogeneous 
solution is preserved. 

The transport equation was derived from the DSE by performing the gradient expansion. This expansion is justified 
only when the quantum and kinetic scales are well separated. We have introduced a mass parameter /i as a separating 
point of these two scales. In ultra-relativistic heavy-ion collisions, the weak coupling condition, g <C 1, is most likely 
to be fulfilled. In this case, low-momentum collective excitations that develop at the momentum scale g\x are well 
separated from typical hard gluon momenta k > \i. We took the difference of the two DSEs, which are in conjugate 
form, and then performed the gradient expansion for the resulting equation under the above conditions. Finally we 
used the gauge covariant GFs, which are obtained by modifying the phase of the conventional GF through Wilson 
links, to derive our final result of the Boltzmann equation (p9|). The sum of the two DSEs and its subsequent gradient 
expansion gives the mass-shell constraint equation for the gauge covariant GF. 

The quantum kinetic equation was shown to be a natural generalization of the classical one, even though that it 
shows a more complicated non-Abelian structure. A notable feature of our quantum result is that, as in the classical 
case, it contains a term that corresponds to the colour precession, the non-Abelian analogue to the Larmor precession 
for particles with magnetic moments in a magnetic field. This term is necessary to guarantee the gauge covariance of 
the quantum kinetic equation. 

The difference between the two conjugate DSEs (Q) and ([|5|) gives the collision part. We have obtained this 
in a gauge covariant form and derived its linearized form in the vicinity of equilibrium. Applying the transversality 
requirement for the gauge covariant Wigner functions, which arises from the background gauge condition in the leading 
order, 0{g 2 T), we have explicitly identified the collision and the damping terms. A similar equation was previously 
derived in Ref. Q in the background Coulomb gauge. However, in the background covariant gauge the results are 
more compact and have an explicit Lorentz covariance. The contribution of the precession term in the kinetic part of 
the Boltzmann equation was shown to be subleading, with respect to the off-equilibrium fluctuations, thus it is not 
there in the vicinity of equilibrium. 

The current approach can be applied to study the propagation of high energy partons (jets) through a hot and cold 
QCD medium. This is because the coherent partons can be treated as the classical background field. Over the past 
few years, substantial progress has been made in understanding the induced gluon radiation in a QCD medium § |. 
However, owing to the complexity of the problem, some idealized and simplified conditions were assumed in order 
to obtain analytical or numerical results. In particular the QCD medium is usually assumed to be in local chemical 
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and thermal equilibrium. We note that following the approach presented in this work, one can study the energy loss 
of the fast parton and the jet quenching via kinetic transport model suitable for computer simulations. In this way 
one could study the influence of the off-equilibrium effects on parton propagation and radiation. We will address this 
problem in the future. 
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APPENDIX A: BARE VERTICES 

In this appendix, we obtain bare vertices by taking derivatives of the classical action with respect to corresponding 
fields: 



S 2 S 



5Qp h {x 1 )5Q^{x 2 ) 
= 6 ht J d i uS i (x 1 - u) g pT dl - d UT d up + ^d UT d up S 4 (x 2 - u) , 

where m — (p, h, x\) and n — (r, i, x 2 ); 



rnnp 



S 3 S 



SQP h (x 1 )SA^{x 2 )SQ^( X3 ) 
= 9f M3 [W& -di) T + g P r (di -d 2 + -d 3 ) 

d 4 uS 4 (xi — u)5 A (x 2 — u)5 i (x 3 — u) , 



+g VT (d 2 -d 3 <9i 



S 3 S 



5QP h (x 1 )SQ ri (x 2 )SQ^(x 3 ) 
= gf hl3 [W& - di) T + g P Adi - d 2 ) v 

+g V r(d 2 - d 3 ) p j d i uS i (x 1 - u)5 i (x 2 - u)S i (x 3 - u) , 
where m — (p, h, xi), n — (r, i, x 2 ), p = (rj,j, x 3 ) and d\ p = dS(xi — u)/du p ; 



^mnpq(A Q ) 



S A S 



mnpq 



5QP a {x 1 )5Q ub (x 2 )5A Xc (x 3 )8A' 7d (x i ) 
- g 1 ! lda ! lch {g a \g pv - g^g P x + -g\»g^ 

~g 2 f lca f ldh (g a \g P v - g P ag v \ + -g^g^x) 
-g 2 I lcd f la \g P xg.u - g P *g„x) 

J d A u5 i {xi — u)8 4 (x 2 — u)8 4 (x 3 — u)(5 4 (x4 — u) 

5^ 

SQ^ a (x 1 )8Q vb (x 2 ) 5Q Xc (x 3 ) 5A ad (x 4 ) 

-g 2 f lda f lcb (g*xg P „-9»„g P x) 
-g 2 f lca .f ldb (g*xg P u - g^g.x) 
-g 2 f lcd f la \g P xg*v - g P «g„x) 

J d i uS i (xi — u)S 4: (x 2 — u)S i (x 3 — u)S i (x4 — u) , 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where m = (p,a,xi), n = (u,b,x 2 ), p = (A, c,x 3 ) and q = (a,d,Xi). We can prove that rSn P9 (Q 4 ) = r^ P9 (AQ 3 ); 



5 2 S 



5C i (x 2 )5Cr(x 1 ) 



J d A u5 i (x\ — u)d 2 8 A (x 2 — u) 



(A6) 



26 



where m — (h, x\) and n = (i, x 2 ); 

S 3 S 



T^UCCQ) 



5C*(x 2 )5C { Xl )SQpi{x 3 ) 



gf h ^ J d 4 u5 4 ( Xl - u)d up 5\x 2 - u)6 4 (x 3 - u) , (A7) 



T^UCCA) 



mnj>\ J h 

SC i (x 2 )6C (xi)5APJ(x 3 ) 



= 9.f hjl J d 4 u{s 4 ( Xl - u)d up [s 4 (x 2 - u)5 4 {x 3 - u) 

+6 4 (x 1 - u)S 4 (x 3 - u)d up 5\x 2 -u)}, (A8) 
where m — (h,x\), n = {i,x 2 ) and p = (j, p,x 3 )] 



TW(CCAQ) = 



6C l (x 2 )SC h { Xi )SApi {x 3 )6Q*>k (x 4 ) 



= 9 2 ! h ° a l akl 9 P n J #u5\x x - u)5 4 (x 2 - u)S 4 (x 3 - u)S 4 (x 4 - u) , (A9) 
where m — (h, x\), n = (i, x 2 ), p = (j, p, x 3 ) and q = (fc, 77, X4). 

r { X q (ccA") = 

5C l {x 2 )5C {x 1 )SApi(x 3 )SA^{x 4 ) 
= 9 2 (.f hla .f akt + .f hka .f ari )9pv 

■ J d 4 uS 4 ( Xl - u)S 4 (x 2 - u)S 4 (x 3 - u)5 4 {x A - u) , (A10) 

with m = (h, xi), n — (i, x 2 ), p = (j, p, x 3 ) and q = (k, 77, X4). 



APPENDIX B: RELATIONS AND IDENTITIES FOR GREEN FUNCTIONS 

Relations between the GF and the CGF for the gluon field Q are given by 

S 2 W 5W 6W 

«W><> = J^ q + JT p JTr G ^ Q) + iQMh 

(Q Q Q ) — + S ' — — + — — + — — 

i 2 SJ m SJ n SJ p iSJ m SJ n SJ p idJ n 5J p 5J m i5J m 5J p dJ n 

= G mnp (Q 3 ) + G mn (Q 2 )(Qp) + G np {Q 2 )(Qm} + G m p(Q 2 )(Qn) , 

(Q Q Q Q) — + - — — 

i 3 5J m 5J n 5J p 8J q i 2 5J m SJ n 8J p SJ q 

S 3 W 5W | 6 3 W 5W 

i 2 5J m dJn8J q SJp i 2 5J m 5J p 5Jp SJ n 

6 3 W SW S 2 W S 2 W 

i 2 SJ n SJpSJ q SJ m iSJ m SJ n idJpdJ q 
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5 2 W S 2 W 5 2 W S 2 W 



iSJ m SJp iSJ n 5J q iSJ m 5J q iSJ n SJp 

G m npq{Q ) "i" G mn p(Q ){Qq) 

+G mnq (Q 3 ){Q p ) + G npq (Q 3 )(Q m ) 
+G mpq (Q 3 )(Q n ) + G mn {Q 2 )G pq (Q 2 ) 

+G mp (Q 2 )G nq (Q 2 ) + G mq (Q 2 )G np (Q 2 ) , (Bl) 

where G mn (Q 2 ), G mnp (Q 3 ) and G mnpq (Q 4 ) are the 2-, 3- and 4-point CGFs for Q. The above equations are easy to 
prove if we recall, for example, that {Q p Q q } = Z~ 1 5 2 Z/(i 2 5J p 5J q ) and W = —i\nZ. For convenience, we sometimes 



T mnp{ 

The relations between the GF and the CGF for Q and C/C are: 



use a short hand notation (mn) = G mn (Q 2 ), (mnp) = G mnp (Q 3 ), etc. 



— S 2 W — 

(C m G n ) — . — — = G m n(CC) , 



{C rnCn Qp) 



S 3 W S 2 W sw 



i 2 S^ n S^ m SJ p i5£ n 8£ m 8J P 

GmnpiCCQ) + G m n(CC)(Q p ) 



CQQ ^ _ ^W | S 3 W SW 

i 3 Sl n S£_ m SJ p SJ q i 2 5£ n S£, m 5J p SJ q 
5 3 W SW S 2 W 5 2 W 



+ 



i 2 5iJi m 5J q SJ P i8£ n 8£ m iSJ p 6J q 

= G m npq(CCQ ) + G m np(CCQ)(Q q ) 

+G mnq (CCQ)(Qp) + G mn (CC)G pq (Q 2 ) , (B2) 

where G mn (CC), G mnp (CCQ) and G mnpq (CCQ 2 ) are the 2-, 3- and 4-point CGFs for Q and C/C. We also use the 
following simplified notation: ([mn]) = G mn (CC), ([mn]p) = G mn p(CCQ) and ([mn}pq) = G mn pq(CCQ 2 ) etc. 
A higher-rank CGF is related to lower-rank ones and 1-PI vertices through the following identities 

(mnp) = iY m > n / p /(Q 3 )(m'm)(n'n)(p'p) , 
(mnpq) = ir m i n , p < q/ (Q A )(m'm)(n'n)(p'p)(q'q) 
+iT m >n' P ' (Q 3 ){m'mq) (n'n) (p'p) 



(Q 3 ) (m'm) (n'nq) (p'p) 



{Q 3 ){m'm){n'n)(jp'pq) , (B3) 



(\mn]p) = iT n i m i p i (CCQ){[mm'])([n 'n])(p 'p) , 
([mn]pq) = iT n > m > p > q >(CCQ 2 )([mm'])([n'n])(p'p)(q'q) 
+iT n > m > P > (CCQ) ( [mm'} q){n'n) {p'p) 

n' m' p' 

(CCQ)([mm'])([n'n}q)(p'p) 
+iT n , m , p ,{CCQ){[mm']){[n'n]){p'pq) , (B4) 

which can easily be derived by repeatedly taking the derivative with respect to (Q) for 5T/5(Q) = — J. 
The 4-point GF can be expressed in terms of the lower-rank CGFs as follows: 

(Qn'Qp'QnQp) = {n'p'){np) + {n'n)(p'p) + (n'p){p'n) + (n'p'np) 
= (n'p')(np) + (n'n)(p'p) + (n'p)(p'n) 
+ir r , s , t , u , (Q 4 ) (r V) (s'p') (t'n) (u'p) 
+iT r , s , t ,{Q 3 ){r'n'p){s'p'){t'n) 
+iT r , s , t ,{Q 3 ){r'n'){s'p'p){t'n) 

+iT r , s , t ,{Q 3 ){r'n'){s'p'){t'np) . (B5) 
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Note that in the second equality the fifth and sixth terms are equal. We can easily identify in (Qn'Qp'QnQp) the 
disconnected Green function (n'p')(np). After dropping it, we get the connected part (Qn'Qp'QnQp) c- 
The 5-point GF can be expressed in terms of the lower-rank CGFs as follows: 

l^l plq ,(Q 4 )(Qn'Q P 'Q q 'QnQp) = \r { l, p , q ,(Q A )[{n'p' q 'np) + (n'p')(q'np) 

+(n q )(p np) + (p q )(n np) + (n n)(p q p) 
+(p n)(n q p) + (q n)(npp) + (n p)(p q n) 

+ (p'p)(n'q'n) + (q'p)(n'p'n) + (np)(n'p'q) , (B6) 

where the last term is a disconnected one. There are three groups of terms that are equal, respectively: (2nd, 3rd, 
4th), (5th, 6th, 7th) and (8th, 9th, 10th). After collecting these terms, we get 

l^l' P ' q '(Q 4 )(Qn'Qp'Q,'QnQp) c = lr^ n>p , g ,(Q 4 )[3(n'p')(q'np) + 3(n'n)(p'q'p) 

+3(n'p)(p'q'n) + (n'p'q'np)] , (B7) 

where the 5-point CGF can be expanded as 

{n'p'q'np) = iT r , s , t , u , v ,(Q 5 )(r'n')(s'p')(t' q ')(u'n)(v'p) 
+iiTr' s ' t 'u' (<3 4 ) (r'n'p) (s'p') (t'q')(u'n) 
+iiTr' 8 ' t 'u' (<3 4 ) (r'n'n) {s'p')(t'q'){u'p) 
+*r rVt , u , (Q 4 )(r'n')(s'p')(t'q') (u'np) 
+6iT r , s , t ,(Q 3 )(r'n'n)(s'p'p)(t' q ') 

+3ir r , s , t ,(Q 4 )(r'n'np)(s'p')(t'q') . (B8) 
The 4-point GF {C m >C n 'Q n Qp) can be expressed in terms of a lower-rank CGFs as follows: 

-^m' n ' m (CC<3) (Cm' C n > Q n Qp) 

= V^ n , m (CCQ){([m'n']np) + ([m'n'])(np)} 

= r^ m (CCQ){^ r , t , u ,(CCQ 2 )([mV])([ S V])(^)(nV) 
+^^ sV , ^ ,(CCO)([TOV>)([ S V])(^'n)+^^ sVt '(CCO)([TO , r'])([ S 'n>)(^'n) 
+iT s , r , t ,(CCQ)([m'r'])([s'n'])(t'np) + ([m'n'])(np)} , (B9) 

where the last term is the disconnected one. 

The 5-point GF {C m >C n 'Qp'QnQp) can be written as 

^rn'n'p'm(CCQ 2 ) (C ' m i C n > Qp' QnQp) 

= T^ p , m (CCQ 2 ){([m'n']n)(p'p) 

+([m'n']p){p'n) + {[m'n']p'){np) + {[m'n']p'np) + {[m'n'])(p'np)} , (BIO) 

where the fourth term is the disconnected one, and the 5-point CGF G m 'n' P 'np(CCQ 3 ) is given by: 

([m'n']p'np) = iT [s , r , ]t , u , v ,([m'r'})([ S 'n'})(t'p')(u'n)(v'p) 
+iT [s , r , ]t , u , {[m'r']p) ([s'n 1 ] ) (t'p')(u'n) + iT [s , r , ]t , u , ([m'r'] )([s'n']p) (t'p') (u'n) 
+iT [s , r , ]t , u , ( [m'r'} ) ( [s'n 1 ] ) (t'p'p) (u'n) + iT [s , r , ]t , u , ( [m'r'] ) ( [s'ri] ) [t'p') (u'np) 
+iT [s , r , ]t , u , (\m'r']n) ([s V] ) (t'p')(u'p) + iT [s , r , ]t , u ,([m'r'] )([s'n']n) (t'p') (u'p) 
+iV [s , r , ]t , u , ( [m'r 1 ] ) ( [s'n'} ) (t'p'n) (u'p) 

+i^[s'r']t'([m'r']np)([s'n'])(t'p') + iT [s/r > ]t >([m'r']n)([s'n']p)(t'p') 
+iT [s , r , ]t ,([m'r']n)([ S 'n'])(t'p'p) + iT [s , r , ]t ,([m' r']p)([ S ' n']n)(t' p') 
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+zr [sV , ]t ,([™V])([ S 'n>p)(tV)+zr [sV , ]t ,(Kr'])([ S 'n»(tVp) 
+ir [s , r , ]t ,([mV / ]p)([ S , n'])(tVn)+ir [sV , ]t ,([m , r'])([ S / n']p)(iVn) 

+ir [sV , ]t ,([m'r'])([ S 'n'])(tVnp) , (Bll) 



where we use the following short-hand notations: 



F[s'r']t'u'v' = ^s'r't'u'v'iCCQ 3 ) , 

r[ s 'r']t'u' = r s ' r 't'«' (ccq ) , 
r[ S ' r ']t' = r s / r / t /(ccQ) . (B12) 
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FIG. 3: Connected part of 4-point GF {Q1Q2Q3Q4) , corresponding to Eq. (B5). 



FIG. 4: Connected part of 5-point GF {Q1Q2Q3Q4Q5), corresponding to Eq. (B7) 




FIG. 5: Connected part of 4-point GF (CCQ1Q2}, corresponding to Eq. (B9). 




FIG. 6: Connected part of 5-point GF {CCQ1Q2Q3), corresponding to Eq. ( BIO ) 



